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We report STAR results on the azimuthal anisotropy parameter V2 for strange particles K%, A 
and A at mid-rapidity in Au+Au collisions at ^/s NN = 130 GeV at RHIC. The value of V2 as a 
function of transverse momentum of the produced particles p t and collision centrality is presented 
for both particles up to pt ~ 3.0 GeV/c. A strong pt dependence in V2 is observed up to 2.0 GeV/c. 
The i>2 measurement is compared with hydrodynamic model calculations. The physics implications 
of the pt integrated V2 magnitude as a function of particle mass are also discussed. 

PACS numbers: 25.75.Wd, 25.75.Ld 



Measurements of azimuthal anisotropics in the trans- 
verse momentum distribution of particles can probe early 
stages of ultra-relativistic heavy-ion collisions § g, §. 
In high-energy nuclear collisions, the initial geometric 
anisotropy is established from the overlap between the 
colliding nuclei. The time necessary to build up this spa- 
tial anisotropy is believed to be short because the col- 
liding nuclei are highly Lorentz contracted in the center- 
of-mass system and pass through each other at approxi- 
mately the speed of light. During a ~ 5-50 fm/c period, 
rescattering transfers the initial spatial anisotropy into 
a momentum anisotropy. This momentum anisotropy 
manifests itself most strongly in the azimuthal distribu- 
tion of transverse momenta. The extent to which the 
initial spatial anisotropy is transformed to the measured 
momentum anisotropy depends on the initial conditions 
and the dynamical evolution of the system. In partic- 
ular, anisotropy measurements for nucleus-nucleus colli- 
sions at RHIC energies may provide information about a 
partonic stage that may exist early in the collision evo- 
lution§|||0,|. 



The transverse momentum distribution of particles can 



be described in the form: 

d 2 N dN .„ . ,.. 

where pt is the transverse momentum of the particle, 
4> is its azimuthal angle with respect to the reaction 
plane j^, [l^] and the harmonic coefficients, v n , are 
anisotropy parameters. The second coefficient v-i is called 
elliptic flow, where flow denotes collective behavior with- 
out necessarily implying a hydrodynamic limit. Recent 
experimental results from RHIC fllj, [jj], [if, |J include 
measurements of V2 as a function of collision centrality 
and pt for charged particles with pt < 2.0 GeV/c, and for 
identified charged pions, kaons and protons for p t up to 
~ 0.8 GeV/c. The degree of the anisotropy transfer from 
position to momentum distribution depends on the den- 
sity of the system during its evolution and the scattering 
cross sections of the particles involved (parton and/or 
hadron). As a result, recent theoretical work attempted 
to deduce the initial gluon density from partonic energy 
loss H, and the equation of state from hydrodynamic 
model calculations f7|. 

Most of the anisotropic flow parameters measured to 
date are for non-strange particles [jll], [l2], |l5|, [l7|, 
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|l8| , fHjL Of the studies for identified strange parti- 
cles |L2|, |0|,!l], H, ||> 11' H most have been at much 
lower collision energies. Moreover, previous measure- 
ments of strange particle flow correspond to directed 
flow, i.e. the coefficient v\. At the CERN SPS, quan- 
titative differences between multi-strange baryons and 
non-strange hadrons were observed in transverse radial 
flow in Pb + Pb collisions at ^/s^ = 17 GeV ge| |7|. 
A physical scenario in which multi-strange baryons do 
not participate in a common expansion and thus decou- 
ple early from the collision system due to their small 
hadronic cross sections, was proposed to explain this ob- 
servation pfi|j . This explanation suggests that it may be 
possible to obtain insight into very early stages of the col- 
lisions by studying the elliptic flow of strange particles. 

In this paper, we report the first measurement of 
the azimuthal anisotropy parameter V2 for the strange 
particles Kg, A and A from Au + Au collisions at 
^/s NN = 130 GeV. Our measurement of V2 for differ- 
ent centralities as a function of pt using the Solenoidal 
Tracker At RHIC (STAR) extends to a p t of about 
3.0 GeV/c, much higher than previously measured for 
identified charged pions, kaons and protons pjl . 

The STAR detector , with its azimuthal symmetry 
and large acceptance, is ideally suited to measure ellip- 
tic flow. The detector consists of several sub-systems in 
a large solenoidal magnet. For collisions in its center, 
the Time Projection Chamber (TPC) measures charged 
tracks in the pseudo-rapidity range \r]\ < 1.5 with 2ir az- 
imuthal coverage. During the year 2000 data taking the 
STAR magnet operated at a 0.25 Tesla field, allowing 
tracking of particles with p t > 0.075 GeV/c. A scintil- 
lator barrel surrounding the TPC, the Central Trigger 
Barrel (CTB), measures the charged particle multiplic- 
ity (for triggering) from within \rj\ < 1. Two zero-degree 
calorimeters J3(j located at ±18.25 m from the nominal 
interaction region, sub-tending an angle 8 < 0.002 radi- 
ans, primarily detect fragmentation neutrons. Two ZDCs 
in coincidence provide a minimum-bias trigger and the 
CTB is used for a central trigger. This analysis uses 
201 x 10 3 minimum-bias and 180 x 10 3 central events. 

We reconstructed both Kg — ► n + + n~ and A (A) — > 
p + 7r~(j? + 7r + ) from their charged daughter tracks de- 
tected in the STAR TPC ||. Using the energy loss of 
the charged tracks in the TPC gas, we select candidates 
for protons, anti-protons and pions. The mass and the 
kinematic properties of the neutral particle candidates 
are extracted from the decay vertex and daughter particle 
kinematics. Fig. |l| shows the invariant mass distributions 
for 7r + 7r~ showing a Kg mass peak and for pir~ showing a 
A mass peak. The dashed lines are fits to the background 
and the signal. We determined that the background is 
dominated by combinatorial counts by rotating all posi- 
tive tracks 180 degrees in the transverse plane and recon- 
structing the K s and A(A) decay vertices. This proce- 
dure destroys all real vertices in the TPC acceptance so 
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FIG. 1: Invariant mass distributions for tt + iy~ showing a K% 
mass peak (left panel) and for pn~ showing a A mass peak 
(right panel). Fitting results are shown as dashed lines in the 
figure. For presentation a greater number of events has been 
used for the A plot. 

that we can describe the combinatorial contribution to 
the invariant mass distributions. The observed masses, 
496 ± 8 MeV/c 2 for tt+tt" and 1116 ± 4 MeV/c 2 for pit, 
are consistent with accepted values |32| an d the widths 
are determined by the momentum resolution of the de- 
tector. The particles used for the V2 analysis are from the 
kinematic region of \y\ < 1.0 and 0.2 < pt < 3.2 GeV/c 
for K s or 0.3 < p t < 3.2 GeV/c for A + A, where y is 
the particle's rapidity. No significant differences in ellip- 
tic flow are observed between A and A, so because of the 
limited statistics, A and A are summed together. 

We choose the requirements for Kg and A(A) daugh- 
ter candidates to maximize statistics. For Kg, we re- 
quire the daughter tracks to have a distance-of-closest- 
approach (dca) to the collision vertex > 1.0 cm. For the 
A(A) reconstruction, we choose pion-like tracks with a 
dca > 1.5 cm and proton-like tracks with a dca > 0.8 cm. 
We use the peak in the invariant mass distribution to 
measure the yield of Kg or A + A particles for different 
values of </> and pt ■ Using the <fi bin center for the value 
of 4>, we evaluate V2 as a function of p t by calculating 
(cos(2(f>)} for different values of p t . Using the yield to 
calculate V2 = (cos(2(f>)} enables us to measure elliptic 
flow for identified particles beyond the pt region where 
the identification of particles via their energy loss in the 
TPC gas fails @. 

The real reaction plane is not known, but the event 
plane, an experimental estimator of the true reaction 
plane, can be calculated from the azimuthal distribution 
of tracks JO). To determine the event plane, we select 
charged particle tracks with at least 15 measured space 
points, 0.1 < p t < 2.0 GeV/c and \t)\ < 1.0. We also 
require the ratio of the number of space points to the ex- 
pected maximum number of space points for each track to 
be greater than 0.52, suppressing split tracks from being 
counted twice. Events are required to have a primary 
vertex within 75 cm longitudinally of the TPC center. 
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FIG. 2: Elliptic flow Vi as a function of pt for (a) .Kg and (b) 
A + A. Circles and filled squares are for central (0-11%) and 
mid-central (11-45%) collisions, respectively. 
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These cuts are similar to those used in Rcf. and our 
analysis is not biased by them. 

To avoid possible auto-correlations, tracks used for the 
Kg or A(A) reconstruction are excluded from the set of 
tracks used to calculate the event plane. Typically this is 
done by measuring the azimuthal angle between a track 
and the event plane calculated from all other qualify- 
ing tracks within the same event. Then the contribu- 
tion to i>2 from that track is calculated. In this analysis, 
where i>2 is not calculated on a particle by particle basis, 
all tracks that might be used for the reconstruction of 
K s or A (A) are excluded from the event plane calcula- 
tion. Only tracks with a dca < 1.0 cm are used in the 
event plane calculation while the K s vertices don't in- 
clude these tracks. In the A + A analysis all proton-like 
tracks are excluded from the event plane calculation. A 
track is considered proton- like if its energy loss (dE/dx) 
is within three standard deviations of that expected for 
protons. 

When the azimuthal anisotropy is evaluated via V2 = 
(cos(2<fi)), the observed V2 must be corrected to account 
for the imperfect event plane resolution |33| . This reso- 
lution is influenced by two factors that depend on cen- 
trality: the strength of the anisotropy signal and the 
number of tracks used for the event plane calculation. 
We estimate the resolution using the method of random 
subevents fl0|| and use the relative multiplicity, as in 
Ref. ]Tl]| , to measure the event centrality. The maxi- 
mum resolution for the Kg and A + A analysis is found 
to be 0.681 ± 0.004 and 0.582 ± 0.007 respectively and 
is reached in the centrality corresponding to 25-35% of 
the measured cross section. The relatively lower resolu- 
tion for the A + A analysis is caused by the exclusion of 
proton-like tracks from the event plane calculation. 

Elliptic flow as a function of transverse momentum 
for central and mid-central collisions calculated from 
201 x 10 3 minimum-bias and 180 x 10 3 central events 
is shown in Fig. |^. The two particles show a similar p t 
dependence in the two centrality intervals. The pt de- 



FIG. 3: Elliptic flow V2 as a function of pt for the strange 
particles Kg (filled circle) and A + A (open squares) from 
minimum-bias Au+Au collisions. For comparison, Vi of 
charged hadrons (open circles) is also shown. The lines are 
from hydrodynamic model calculations JH|. 

pendence is stronger in more peripheral collisions than 
in the central collisions. A similar dependence was ob- 
served for charged particles in Au + Au collisions at the 
same RHIC energy (l|. 

For this analysis, three main sources contribute to sys- 
tematic errors in the measured anisotropy parameters: 
particle identification, background subtraction, and cor- 
relations unrelated to the reaction plane (non-flow) such 
as resonance decays, jets or Coulomb and Bose-Einstein 
correlations |34|, |3f| . The contribution from the first two 
sources is estimated by examining the variation in vi af- 
ter changing several track and event cuts. We estimate 
that these effects contribute an error of less than ± 0.005 
to V2- The contribution to V2 from non-flow effects, how- 
ever, could be significant, especially in peripheral colli- 
sions. A previous study used the correlation of event 
plane angles from subevents to estimate the magnitude 
of these contributions |36| . Non-flow effects are assumed 
to contribute to the first and second harmonic correla- 
tions by similar amounts, so the magnitude of the first 
harmonic correlation sets a limit on the non-flow contri- 
butions to V2 ■ That study showed that the non-flow sys- 
tematic errors for charged particles are typically +0 and 
-0.005, but are significantly larger in the more periph- 
eral events where the error increases to +0 and -0.035 
for the 58-85% most central events. These estimates are 
confirmed by measurements of V2 using the 4th-order cu- 
mulant method, a method that is insensitive to non-flow 
effects but which leads to larger statistical errors [p7j . We 
assume the systematic errors on V2 for the neutral strange 
particles K s and A + A are similar to those found in the 
analysis of charged particles (l^] . 

To make a comparison with available hydrodynamic 
model calculations, we plot i>2(pt) for both K s and A + A 
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FIG. 4: Integrated elliptic flow V2 as a function of particle 
mass. The gray-band and central line indicates the hydrody- 
namic model results j5j. 



from 201 x 10 3 minimum-bias collisions in Fig. |^. The 
dashed lines represent the hydrodynamic model calcu- 
lations H for (from top to bottom) pions, kaons, pro- 
tons, and lambdas. Also shown in the hgure is V2(pt) 
for charged hadrons |Q. Within statistical uncertainty, 
the Kg results are in agreement with the v 2 of charged 



kaons (not shown) 12 . We observe that v 2 for both 
strange particles increases as a function of pt up to about 
1.5 GeV/c, similar to the hydrodynamic model predic- 
tion. In the higher p t region however (p t > 2 GeV/c), 
the values of v 2 seem to be saturated. It has been sug- 
gested that the shape and height of v 2 above 2-3 GeV/c 
in a pQCD model is related to energy loss in an early, 
high-parton-density, stage of the evolution [|| . 

The pt integrated anisotropy parameters for charged 
hadrons, Kg, and A + A from minimum-bias collisions 
are shown in Fig. ||. The integrated values of v 2 are 
calculated by parameterizing the yield with the inverse 
slope parameter of exponential fits to the K s or A + A 
transverse mass distributions (3^, |3^] . The integrated v 2 
is dominated by the region near the particle's mean p t 
and is insensitive to the upper and lower bounds of the 
integration. Although the V2(pt) of A + A is below the 
V2(pt) of K s for most p t , as shown in Fig. [|, the p t in- 
tegrated v 2 values increase with the particle mass. This 
increase is partly due to the relatively higher mean pt of 
the A + A compared to the Kg. In hydrodynamic models, 
although the spatial geometry of the pressure gradient 
and the resultant collective velocity are the same for all 
particles, massive particles tend to gain larger transverse 
momenta and so develop a larger elliptic flow. The hy- 
drodynamic model calculations (j| , shown as a gray-band 
and central line, are, within errors, in agreement with this 
result. The width of the gray-band in Fig. [l] indicates 



the uncertainties of the model calculation, mostly due to 
the choice of the freeze-out conditions. The increase of 
v 2 with particle mass indicates that significant collective 
motion, perhaps established early in the collision, is an 
effective means to transfer geometrical anisotropy to mo- 
mentum anisotropy. The nature of the particles during 
this process, however, whether parton or hadron, and the 
degree of thermalization for strange particles during the 
collective expansion remains an open issue. 



In summary, we have reported the first measurement 
of the anisotropy parameter, v 2 , for Kg and A + A, from 
Au + Au collisions at ^/s NN = 130 GeV. The v 2 values 
as a function of pt from mid-central collisions are higher 
at each j>i than v 2 from central collisions. Hydrodynamic 
model calculations seem to adequately describe elliptic 
flow of the strange particles up to a pt of 2 GeV/c. For 
Pt above 2 GeV/c, however, the observed v 2 seems to 
saturate whereas hydrodynamic models predict a contin- 
ued increase with p t . The p t integrated v 2 as a function 
of particle mass is consistent with a hydrodynamic pic- 
ture where collective motion, established by a pressure 
gradient, transfers geometrical anisotropy to momentum 
anisotropy. Although the hadronic scattering cross sec- 
tions of strange and non-strange particles may be differ- 
ent, we have yet to see deviations in the measured v 2 
from hydrodynamic calculations at low pt for strange or 
non-strange particles. In a possible partonic phase prior 
to the hadronic epoch, the hadronic scattering cross sec- 
tions for the final hadrons are not relevant. As such, if the 
elliptic flow of identified particles proves to be indepen- 
dent of their relative hadronic cross sections, it may be 
evidence that v 2 is established during a partonic phase. 
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